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Abstract 

We show that transition to longitudinal instability of strong detonation solu- 
tions of reactive compressible Navier-Stokes equations is generically associated 
with Hopf bifurcation to nearby time-periodic "galloping" , or "pulsating" , so- 
lutions, in agreement with physical and numerical observation. The analysis is 
by pointwise semigroup techniques introduced by the authors and collaborators 
in previous works. 
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1 Introduction 

Motivated by physical and numerical observations of time-oscillatory "galloping" or 
"pulsating" instabilities of detonation waves [MT, BMR, FW, MT, AIT, AT, Fl, F2, 
KS], we study Hopf bifurcation of viscous detonation waves, or traveling- wave so- 
lutions of the reactive compressible Navier-Stokes equations. This extends a larger 
program begun in [Zl, LyZl, LyZ2, JLW, LRTZ] toward the dynamical study of vis- 
cous combustion waves using Evans function/inverse Laplace transform techniques 
introduced in the context of viscous shock waves [GZ, ZH, ZS, Zl, MaZ3], continuing 
the line of investigation initiated in [TZl, TZ2, SS, TZ3] on bifurcation/transition 
to instability. 

It has long been observed that transition to instability of detonation waves occurs 
in certain predictable ways, with the archetypal behavior in the case of longitudinal. 
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or one-dimensional instability being transition from a steady planar progressing wave 
U {x, t) = U{xi — st) to a galloping, or time-periodic planar progressing wave U{xi — 
st,t), where U is periodic in the second coordinate, and in the case of transverse, or 
multi-dimensional instability, transition to more complicated "spinning" or "cellular 
behavior"; see [KS, TZl, TZ2], and references therein. 

The purpose of this paper is, restricting to the one-dimensional case, to establish 
this principle rigorously, arguing from first principles from the physical equations 
that transition to longitudinal instability of detonation waves is generically associ- 
ated with Hopf bifurcation to time-periodic galloping solutions. 

1.1 The reacting Navier-Stokes equations 

The single-spccics reactive compressible Navier— Stokes equations, in Lagrangian co- 
ordinates, appear as [Ch] 

dtT - dxU = 0, 
dtu + dxP = dxivT'^dxu), 
dtE + dx{pu) = 5^ [qdr^'^d^z + nr'^d^T -\- uT'^ud^u) , 

{dtz + k4>{T)z = d^{dT-''d^z), 



(1.1) 



where r > denotes specific volume, u velocity, E > total specific energy, and 
< z < 1 mass fraction of the reactant. 
The variable 

U := (r, u, E, z) G 

depend on time t G 1R+, position x G M, and parameters v, k, d, k, q, where i/ > is a 
viscosity coefficient, k> and d> are respectively coefficients of heat conduction 
and species diffusion, k > represents the rate of the reaction, and q is the heat 
release parameter, with q > corresponding to an exothermic reaction and q < 
to an endothermic reaction. 

In (1.1), T = T{t, e,z)>0 represents temperature, p = p{t, e, z) pressure, where 
the internal energy e > is defined through the relation 

1 2 

E = e + -u"^ -\- qz. 

In (1.1), we assume a simple one-step, one-reactant, one-product reaction 

a'^^B, z:=[A], [A] + [B] = 1. 

where (j) is an ignition function. More realistic reaction models are described in 
[GS2]. 

In the variable U, after the shift 

X ^ X — st, s G M, 
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the system (1.1) takes the form of a system of differential equations 

(1.2) dtU + d^{F{U)) = d^{B{U)d^U) + G{U), 

I \ 




V -k<t>{T)z I 



KT-'drT -KUT-'deT + vT-'u KT'^deT kt'^ {d;,T - qdeT) + qdr''^ 
V J 

The characteristic speeds of the first-order part of (1.1), i.e., the eigenvalues of 
duF{U), are 

(1.3) { —5 — a, —s, -s + cr , ^^-s^ }, 

fluid eigenvalues reactive eigenvalue 

where a, the sound speed of the gas, is 

1.2 Assumptions 

We make the following assumptions: 

Assumption 1.1. We assume a reaction-independent ideal gas equation of state, 

p = Tr-^e, T = c~^e, 

where c> is the specific heat constant and T is the Gruneisen constant. 

Assumption 1.2. The ignition function (p is smooth; it vanishes identically for 
T <Ti, and is strictly positive for T > Ti. 

1.3 Coordinatizations 

We let 

w:= {u,E,z) eM.^, v:= {t,u,E) eM.^. 
Then we have the coordinatizations 

U = (V,Z) = iT,w). 

4 



where 



and 



P 

pu 

V / 



VT 



s{e)U, G 








In particular, Assumption 1.1 implies that in the (r, w) coordinatization, B takes 
the block-diagonal form 


b 



B 



where b is full rank for all values of the parameters and U; the system (1.2) in (r, w) 
coordinates is 

{dtT - sdxT - JdxW = 0, 
dfW + dxf{T, w) = dx{b{T, w)dxw) + g{w), 

with the notation 

(1.4) J := ( 1 ) , f:=(pu\-sw, g := ( 

V y V -kHT)z 

In the {v,z) coordinatization, the system (1.2) takes the form 

f dtv + dj\v, z) = d^{b\{v)d,v + b\{v)d^z) 
y dtz — sdxZ + k(p{T)z = dx{dT~'^dxz). 

where the flux is /t* = {—u — ST,p — su,pu — sE), and, under Assumption 1.1, the 
diflFusion matrices are 

b{= I i/r-i , bl= I 

\0 T-^{u - Kc-^)u KT-^c-^ J \ g(dr-2 - Acr-i) 

Note that, in the {v, z) coordinatization, the first component is a conservative 
variable, in the sense that dfV is a perfect derivative, hence 

(1.5) {v{x,t) -v{x,0))dx = 0, 
for v{t)-v{0) e W^^^' 



1.4 Strong detonations 

We prove in this article stability and bifurcation results for viscous strong detona- 
tions of (1.1), defined as follows: 

Definition 1.3. A one-parameter, right-going family of viscous strong detonations 
is a family {U^}^^^ of smooth stationary solutions of (1.2), associated with speeds 
s(e), s{e) > 0, model parameters {v, K,d,k,q){s) and ignition function cj)^, with 
U^,(f)^, {sjU, K,d,k,q){e) depending smoothly on e in L°° x L°° x R^, satisfying 

(1.6) U%x,t) = U%x), lim C7=(x) = C/^, 
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connecting a burned state on the left to an unbumed state on the right, 

(1.7) zL=0, 4 = 1, 

with a temperature on the burned side above ignition temperature 

(1.8) Tf > Ti, 
and satisfying the Lax characteristic conditions 

(1.9) (7- := a{Ui) >s>a+:= a{Ul), 
uniformly in e. 

Consider a standing wave (1.6), U = (r, u, E, z), solution of (1.2), with endstates 
U± = {t±,u±, E±, z±). It satisfies the hnear constraint 

— s(r — r_) = u — 

the system of ordinary differential equations 

UT~^u' = p — SU — {p — SU)-, 

kt~^c~^E' + T~^{y — Kc~^)uu' = pu — sE — {pu — sE)- 

(1.10) { + (kt-^c-i - dT-^)qy, 

z' = y, 

dr^'^y' = -sy + k(j){T)z, 
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and the Rankine-Hugoniot relations 



(1.11) 



-s(r+ — r_) = u+ — U-, 

{p-su)+ = {p-su}-, 
{jpu — sE)^ = (jpu — sE). 

y± = o, 

(l>{T±)z± = 0, 



expressing the fact that {u±,E±,Q,z±) are rest points of (1.10). 

Prom (1.9) and (1.11), we note that the right endstate of a strong detonation 
satisfies 



(1.12) ct>{Tl) = 0, 
which, by Assumption 1.2, imphes also 

(1.13) cP'iTl) = 0. 

Lemma 1.4. Under Assumptions 1.1, 1.2, if q > and s is large enough with 
respect to q, then for any z+ € (0,1], there exists an open subset O- in M^, such 
that any left endstate U- = (v-,0) with v- G 0_ satisfies (1.8) and (1.9), and is 
associated with a right endstate [/+ = (-?;+, z+) satisfying r+ < Tj, (1.9) and (1.11). 

The existence of strong detonations was proved by Gasscr and Szmolyan [GSl] 
for small dissipation coefficients v, k and d. Wc restrict throughout the article to 
strong detonations with left cndstatcs as in the above Lemma. 

Remark 1.5. In the small-heat-release limit q ^ 0, the equations in {y,z) (in sys- 
tem (1.10)) are decoupled from the fluid equations; in particular, strong detonations 
converge to ordinary nonreacting gas-dynamical shocks of standard Lax type, the 
existence of which was first proved by Gilbarg [G] . 

A consequence of Lemma 1.4 is that strong detonations converge exponentially 
to their endstates, a key fact of the subsequent stability and bifurcation analysis. 

Corollary 1.6. Under Assumptions 1.1, 1.2, let {U^}e be a family of viscous strong 
detonations. There exist C,r]o > 0, such that, for k >0 and j G {0, 1}, 

(1.14) ' ' , _ ^' - 
|^a^([/^-[/^)(a;)| <Ce-^°l^l, x > 0. 

In particular, |(C/^)'(x)| < Ce-^ol^l, for all x. 

Remark 1.7. In the ZND limit, strong detonations are transverse orbits of (1.10), 
a result proved in Section 3.6 of [LyZ2], after [GSl]. 

Lemma 1.4 and Corollary 1.6 are proved in Section 2. 
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1.5 Structure of the equations and the profiles 

System (1.2), seen as a system in r, u;, satisfies 

(Al) the convection terms in the equation in r are linear in (t,w); 

(A2) the diffusion matrix b is positive definite. 

For strong detonation waves, the convection terms in (1.1) satisfy 

(HI) The convection coefficient s{e) in the evolution equation in r is nonzero, uni- 
formly in e. 

(H2) The spectrum of duF, given in (1.3), is real, simple, and nonzero, uniformly 
in e. 

System (1.2) satisfies the Kawashima dissipativity condition 
(H3) For ah e, for ah ^ G R, 

at the endstates of a family of strong detonations. In (H3), a denotes spectrum 
of a matrix, and > is independent of and e. To verify (H3), it suffices, by a 
classical result of [ShK], to check that (1.2) has a symmetrizable hyperbolic-parabolic 
structure, and that the genuine coupling condition holds. These conditions are 
coordinates-independent, and easily checked in (r, u, e) coordinates. 
Finally, the assumption 

(H4) Considered as connecting orbits of (1.10), lie in a smooth one-dimensional 
manifold of solutions of (1.10), obtained as a transversal intersection of the 
unstable manifold at Ut and the stable manifold at U^, 

holds in the ZND limit, as stated in Remark 1.7. Under (H4), in a vicinity of U^, the 
set of stationary solutions of (1.2) with limits at ±oo is a smooth one-dimensional 
manifold, given by {ij^{- — c), c G M}, and the associated speed e — s{e) is smooth. 

Conditions (A1)-(A2),(H0)-(H4) are the assumptions of [TZ3] (where G = 0), 
themselves a strengthened version of the assumptions of [MaZ3] . 

1.6 The Evans function 

A central object in the study of stability of traveling waves is the Evans function 
D{£, •) (precisely defined in Section 3.1.7), a Wronskian of solutions of the eigenvalue 
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equation {L{s) — X)U = decaying at plus or minus spatial infinity [AGJ]^, where 
the linearized operator L is defined as 

(1.15) L{e) := -d,{A •) + d,{B{W)d^ ■) + %G(C7^), 
with the notation 

(1.16) A := -duF{U') + {duBiU') -XU')'. 
Recall the important result of [LyZ2]: 

Proposition 1.8 ([LyZ2], Theorem 4). Under Assumptions 1.1 and 1.2, let {U'^}e 
be a one-parameter family of viscous strong detonation waves satisfying (H4) . 

For all e, the associated Evans function has a zero of multiplicity one at X = : 

D{e,0)=0, and D'{e,0)^0. 

Proof. By translational invariance, D{e, 0) = 0, for all e. Generalizing similar results 
known for shock waves [GZ, ZS], there was established in [Zl, LyZl, LyZ2] the 
fundamental relation 

(1.17) D'{e,0)=jS. 

In (1.17), 7 is a coefficient given as a Wronskian of solutions of the linearized 
traveling-wave ODE about U; transversality corresponds to 7 0. In (1.17), 6 
is the Lopatinski determinant 

5 := det ( ( t+ — r_ u+ — u- — E^ + q )*^ ), 

(where rj denote the eigenvectors of duF(UV) associated with outgoing eigenvalues, 
F as in (1.2), and tr denotes transverse matrix^) determining hyperbolic stability 
of the Chapman-Jouget (square wave) approximation modeling the detonation as 
a shock discontinuity. Hyperbolic stability corresponds to 5 7^ 0. See [Zl, LyZl, 
LyZ2, JLW] for further discussion. By (H4), 7 / 0, while (5 / by direct calculation 
comparing to the nonreactive (shock-wave) case. □ 

Remark 1.9. The vectors r^, and r^ correspond to outgoing modes to the left 
of X = 0, see Section 3.1.2 and Figure 4. (The fluid modes rJ, 1 < i < 3, are 
ordered as usual by increasing characteristic speeds: —s — a- <—s<Q<—s + a-, 
so that r^ is incoming.) 

^For applications of the Evans function to stability of viscous shock and detonation waves, see, 
e.g., [AGJ, GZ, ZS, Zl, LyZl, LyZ2, LRTZ]. 

^This notation will be used throughout the article. 
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1.7 Results 



Let X and Y be two Banach spaces, and consider a traveling wave U solution of a 
general evolution equation. 

Definition 1.10. A traveling wave U is said to be X linearly orbitally stable 

if, for any solution U of the linearized equations about U with initial data in X, 
there exists a phase shift 6, such that \U{-,t) — S{t)U'{-)\Y is bounded for <t < oo. 

It is said to be X ^ Y linearly asymptotically orbitally stable if it is X Y 
linearly orbitally stable and if moreover \U{-,t) — 5{t)U'{-)\Y — > as t — > oo. 

Definition 1.11. A traveling wave U is said to be X ^ Y nonlinearly orbitally 
stable if, for each S > 0, for any solution U of the nonlinear equations with 
|f/(-,0) — U\x sufficiently small, there exists a phase shift 5, such that \U{-,t) — 
U{- -6{t),t)\Y <SforO<t<oo. 

It is said to be X ^ Y nonlinearly asymptotically orbitally stable if it is X ^ Y 
nonlinearly orbitally stable and if moreover \U{-,t) — U{- — S(t),t)\Y — > as t ^ oo. 

1.7.1 Stability 

Our first result, generalizing that of [LRTZ] in the artificial viscosity case, is a char- 
acterization of linearized stability and a sufficient condition for nonlinear stability, 
in terms of an Evans function condition. 

Theorem 1.12. Under Assumptions 1.1, 1.2, let {U^}e be a one-parameter family 
of viscous strong detonation waves. 

For all e, is ^LP ^ LP linearly orbitally stable if and only if, for all e, 

(1.18) the only zero of D{£, •) in 5ftA > is a simple zero at the origin. 

If (1.18) holds, tj^ is n H'^ ^ fl linearly and nonlinearly orbitally stable, 
and H n asymptotically orbitally stable, for p > 1, with 

(1.19) \U'{;t) - W{- - 5mLp < CPS - U'\L^nHs{l + tr-^^'-'^\ 
where is the solution of (1.2) issued from Uq, for some 6{-) satisfying 

\S{t)\<C\Ul-U'\LinHS, 
m\<C\Ul-U'\LinHs{l + trl 

1.7.2 Transition from stability to instability 

Tlieorem 1.13. Under Assumptions 1.1, 1.2, let {U^}^ be a one-parameter family 
of viscous strong detonation waves satisfying (H4). 
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Assume that the family of equations (1.2) and profiles undergoes transition 
to instability at e = in the sense that is linearly stable for e < and linearly 
unstable for e > 0. 

Then, one or more pair of nonzero complex conjugate eigenvalues of L{£) move 
from the stable (negative real part) to the neutral or unstable (nonnegative real part) 
half-plane as e passes from negative to positive through £ = 0, while A = remains 

a simple root of D{£, ■) for all e. 

That is, transition to instability is associated with a Hopf-type bifurcation in 
the spectral configuration of the linearized operator about the wave. 

Proof of Theorem 1.13. By Theorem 1.12, transition from stability to instability 
must occur through the passage of a root of the Evans function from the stable 
half-plane to the neutral or unstable half-plane. However, Proposition 1.8 implies 
that D has a zero of multiplicity one at the origin, for all e, and so no root can pass 
through the origin. It follows that transition to instability, if it occurs, must occur 
through the passage of one or more nonzero complex conjugate pairs A = 7 ± ir, 
r 7^ 0, from the stable half-plane (7 < for £ < 0) to the neutral or unstable 
half-plane (7 > for £ > 0). □ 

Our third result and the main object of this paper is to establish, under appro- 
priate nondegeneracy conditions, that the spectral Hopf bifurcation configuration 
described in Theorem 1.13 is realized at the nonlinear level as a genuine bifurcation 
to time-periodic solutions. 

Given A; G N and a weight function a; > 0, define the Sobolev space and associated 
norm 

(1.20) Ht := {/ G S'm, Jf e H\R)}, ||/||^. := ||a;^/||^.. 

Let oJ G be a growing weight function such that, for some ^0 > 0, C > 0, for 
all x,y, 



(1.21) 



1 < u;(x) < e^"(l+l"l')^ 
\u}'{x)\ + \co"{x)\ < Cuj{x), 

u!{x) < Cbj{x — y)oj{y). 



Theorem 1.14. Under Assumptions 1.1, 1.2, let {U^}e be a family of viscous strong 
detonation waves satisfying (H4). 

Assume that the family of equations (1.2) and profiles undergoes transition 
from linear stability to linear instability at e = 0. 

Moreover, assume that this transition is associated with passage of a single com- 
plex conjugate pair of eigenvalues ofL{e), A-i-(£) = 7(£)-|-ir(£) through the imaginary 
axis, satisfying 

(1.22) 7(0) = 0, t{0)^0, d7/d£(0)/0. 
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Then, given a growing weight satisfying (1.21) with 9q sufficiently small, for 
r > sufficiently small and C > sufficiently large, there are functions r ^{f), 
r — > T(r), with e(0) = 0, r(0) = 27r/r(0), and a family of time-periodic solutions 
U^{x,t) of (1.2) with £ = e{r), of period T (r) , with 

(1.23) C-^r < \\U'' - U'\\h2 < Cr, 

Up to translation in x, t, these are the only time-periodic solutions nearby in \\ ■ ||^2 
with period T € [Tq, Ti] for any fixed < Tq < Ti < +oo. 

That is, transition to linear instability of viscous strong detonation waves is 
"generically" (in the sense of (1.22)) associated with Hopf bifurcation to time- 
periodic galloping solutions, as asserted in the title of this paper. 

The choices u = 1 and u = e^^'^^+l^l^^^ are allowed in (1.21), as well as w = 
(1 + |a;p)P, for any real p > 0. In Theorem 1.14, we need, in particular, < 
where r/o is as in Corollary 1.6, so that the spatial localization given by (1.23) is 
less precise than the spatial localization of the background profile U^. The smallness 
condition on is described in Remark 5.9. 

1.7.3 Nonlinear instability 

Wc complete our discussion with the following straightforward result verifying that 
the exchange of linear stability described in Theorem 1.14, as expected, corresponds 
to an exchange of nonlinear stability as well, the new assertion being nonlinear 
instability for £ > 0. 

Theorem 1.15. Under the assumptions of Theorem 1.14, the viscous strong det- 
onation waves undergo a transition at e = from nonlinear orbital stability to 
instability; that is, is nonlinearly orbitally stable for e < and unstable for e > 0. 

1.8 Discussion and open problems 

This analysis in large part concludes the one-dimensional program set out in [TZ2]. 
However, a very interesting remaining open problem is to determine linearized and 
nonlinear stability of the bifurcating time-periodic solutions, in the spirit of Section 
4.3. For a treatment in the shock wave case with semilinear viscosity, see [BeSZ]. 
Likewise, it would be very interesting to carry out a numerical investigation of the 
spectrum of the linearized operator about detonation waves with varying physical 
parameters, as done in [LS, KS] in the inviscid ZND setting, but using the viscous 
methods of [Brl, Br2, BrZ, BDG, HuZ] to treat the full reacting Navier-Stokes 
equations, in order to determine the physical bifurcation boundaries. 

Other interesting open problems are the extension to multi-dimensional (spin- 
ning or cellular) bifurcations, as carried out for artificial viscosity systems in [TZ2], 
and to the case of weak detonations (analogous to the case of undercompressive 
viscous shocks; see [HZ, RZ, LRTZ]). 
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The strong detonation structure considerably simplifies both stability and bifur- 
cation arguments over what was done in [LRTZ]. Wc remark that, at the expense of 
further complication, nonlinear stability of general (time-independent) combustion 
waves, including also weak detonations and strong or weak deflagrations, may be 
treated by a combination of the pointwise arguments of [LRTZ] and [RZ]. 

We remark finally that the restriction to a scalar reaction variable is for simplicity 
only. Indeed, the results of this article (as well as the results of the article by Lyng 
and Zumbrun [LyZ2] from which it draws) are independent from the dimension of 
the reactive equation, so long as the reaction satisfies an assumption of exponential 
decay of space-independent states (with temperature at — oo above the ignition 
temperature). 

Plan of the paper. Lemma 1.4 and Corollary 1.6 are proved in Section 2. 
We give a detailed description of the low-frequency behavior of the resolvent kernel 
for the linearized equations in Section 3, following [MaZS]. In Section 4, we prove 
Theorem 1.12, while Section 5 is devoted to the proof of Theorem 1.14. Finally, in 
Section 6, we prove Theorem 1.15. 

Acknowledgement. Thanks to Bjorn Sandstede and Arnd School for their 
interest in this work and for stimulating discussions on spatial dynamics and bi- 
furcation in the absence of a spectral gap. Thanks to Gregory Lyng for pointing 
out reference [Ch]. B.T. thanks Indiana University for their hospitaUty during the 
collaborative visit in which the analysis was carried out. B.T. and K.Z. separately 
thank the Ecole Polytechnique Federale de Lausanne for their hospitality during 
two visits in which a substantial part of the analysis was carried out. 



2 Strong detonations 

Proof of Lemma 1.4. Let U- be a given left endstate, with = 0, satisfying (1.8) 
and (1.9). We look for a right endstate ?7+, with 2:+ G (0, 1], that satisfies (1.11), 
(1.9), and r+ < Tj. We note that (l.ll)(i) determines Uj^ and that T+ < Tj entails 
(l.ll)(v). 

The Rankine-Hugoniot relations in the (r+,p+) plane are 

'p = —s^T + ci (R), 
p = (co - sr(l + r-^))-^(c2 + sqz+ + ^sV^ - s%r) (H), 

where (R) is the Rayleigh line, corresponding to (1.11) (ii), (H) the Hugoniot curve, 
corresponding to (1.11) (iii), and where 

cq:=U- + st-, ci := P- + s^T-, C2 := {p-U- - sE^) + l-c^s 
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depend on parameters C/_ and s. The temperature and Lax constraints for both 
enstates are 

r < cTTi < T.p. (T)i, 

\T+ip+<(r + i)-is2<rrV (L)^. 

We restrict to left endstates satisfying in the large s regime 

(2.1) r_ = 0(1), p_ = 2s2r-V_ u_ = SU-, 

with ■£(_ = 0(1) andp- = 0{1). Under (2.1), conditions (T)_ and (L)_ are satisfied 
as soon as s is large enough. The Hugoniot curve takes the form 

PH = {u- + r_ - (1 + r-^)Ty\^s^iT - (1 + 2r-i)T_)(r - (1 - 2r-i)r_ - 2u_) 

+ sqz+j . 

Assume that n_ is such that 

(2.2) < (1 - 2r-i)r_ + 2u_ < (1 + 2r-i)r_. 

For any such u_, any given r_ and any g > 0, if s is large enough then, for any 
2+ G (0, 1] , the Hugoniot curve has two zeros r < r , with asymptotic expansions 

(2.3) r = (1 - 2r-^)r_ + 2tx_ + 0{s-^). 

(2.4) r = (1 + 2r-^)r_ - ,-2 + + o(,-3). 

If s is large, by (2.2), tq < t < r, where tq := cqs~^{1 + r~"'^)~^ is the pole of (H). 

The Rayleigh line and the Hugoniot curve have at least one intersection point 
to the right of tq if 

Pr{t) < < pr(t). 

Under (2.2), the inequality < pr{t) holds, and pr{t) < holds as well if in 
addition 

(2.5) P-<- 



2r-i' 



U- 



Let r+ be an intersection point of (R) and (H) to the right of tq. Condition (T)_ 
is satisfied if 

(2.6) r+ = (1 + 2r-i)r_ + s-^f+ + Ois'^), 
with 

(2.7) (1 + 2r^i)r_(p_ - f+) < cTTi. 
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Condition (L)^ is satisfied if 

(2.8) (l + 2r-i)r_ < (l + (l + r)-i)r+, 

which holds under (2.6), if s is large. We plug the ansatz (2.6) in the equation 
PH = PR, to find 

(2.9) f+ = - — E^__((i + r-i)(i + 2r-i)-i) + — ^^^^i— . 
V J + (i + 2r-i)T_^^ ' ' (i + 2r-i)r_ 

The intersection point r+ is an admissible right specific volume if ph(7"+) > and 
PR(r+) > 0. These inequalities holds if 

(2.10) r <r+ < (Q; + l)r_ + s~^p_. 

The inequalities (2.5), (2.7) and (2.10) are constraints on T_,p_, and u_. The lower 
bound on r+ in (2.10) is satisfied in the regime (2.1) if s is large. If we let 

r_ 

then (2.5) holds. Finally, if r_ satisfies 

i< (3 + 2r-i - (1 + 2r-i)r_) < 1 + 

4r_ ^ ' qz+ 

then the upper bound in (2.10) and (2.7) hold as well. 

The Rayleigh line (R), the Hugoniot curve (R) and the temperature (T) and Lax 
(L) constraints are pictured on Figure 2. The black dots represent the intersection 
points of (R) and (H). Note that (L) and (R) imply r_ < r+ for a strong detonation, 
so that only the intersection point to the right to r_ is admissible. (The other 
intersection point corresponds to a deflagration, see for instance [LyZ2], Section 
1.4.) □ 

Proof of Corollary 1.6. Rewrite (1.10) as U' = ^{e,U). Let be the endstates of 
a family of strong detonations. The linearized equations at are governed by 
matrices 



(2.11) dudie, U-_ 







The block triangular structure is a consequence of Assumption 1.1, (1.7), (1.12), 
and (1.13). Under Assumption 1.1, the eigenvalues A of a£ 

^/ _ f VT~'^ \ / duP - S - S~'^drP deP \ 

^' Y T^^{y — K,c'^)u Kc'^T~^ J \ u{duP — s~^dTp) + p udeP — s y ' 
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satisfy 

(2.12) + {sKC-^T^^ + s-Vr±Hs2 + {drp)±))\ + Kc-Vr±2(.s2 - 4) = 0. 

The Lax condition (1.9) implies that the center subspace on both sides is trivial, 

f f 
that the eigenvalues of a_ have opposite signs, and that the eigenvalues of are 

negative. The eigenvalues A of 

r _ ( -sd-^ kd-^(j){T±) \ 

""^'-y 1 ; 

satisfy 

dX^ + s\- k(t){T±) = 0. 

They are non zero and have distinct signs on the — oo side. On the +oo side, there 
is one negative eigenvalue, and a one-dimensional kernel. In particular, is a 
hyperbolic rest point of the linearized traveling-wave ordinary differential equation, 
which implies (1.14) (i) with j = 0, by standard ODE estimates. However, the 
linearized traveling-wave equations at C/_|_ have a one-dimensional center subspace, 
which a priori precludes exponential decay (1.14). 

Prom Lemma 1.4, if U- G O, then the system (1.10) has a line of equilibria that 
goes through [7^. Any center manifold of (1.10) at f/^ contains all equilibria, so by 
dimension count it must consist of equilibria. Therefore, the 4-dimensional stable 
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center manifold at consists (again by dimension count) of the union of the stable 
manifolds of all equilibria. Since solutions off of stable center manifold do not stay 
for all time in small vicinity of center manifold, any traveling-wave orbit must lie 
on the center-stable manifold, so lies on the stable manifold of some equilibrium. 
Exponential decay, (1.14) (ii), j = 0, now follows by the stable manifold theorem. 

To prove (1.14) with j = 1, consider now the traveling-wave ODE in {U,d^U). 
The rest points satisfy 

(2.13) ;?(£,[/) = 0, d,'S{e,U) + du^{e,U)deU = {). 

The kernel of du^is, U^) being one-dimensional, (2.13) has a two-dimensional man- 
ifold of solutions. Let {U^,V_f_) be such a rest point. The linearized equations at 
{U^,V^) are governed by matrices 

f dud{e,Ul) \ 

where the bottom left entry depends on second derivatives of ^. In particular, the 
linearized equations have a two-dimensional center subspace. We can thus argue 
as above that any center manifold consists entirely of equilibria, and that (1.14) (ii) 
holds with j = I. The proof of (1.14) (i) with j = 1 is similar. □ 

3 Resolvent kernel and Green function bounds 

The linearized equations about a traveling wave solution of (1.2) are 

(3.1) dtU = L{e)U, 

where L{e) is defined in (1.15). The coefficients of L{e) arc asymptotically constant 
at zbcxD. Let L±{£) be the associated constant-coefficient, limiting operators: 

L±(£) :=-A±d^ + B±dl + G±, 

with the notation A± := duF{U^), B± := B{Ul), G± := duG{Ul). Let L{e)* de- 
note the dual operator of L{£). Its associated constant-coefficient, limiting operators 
are L±{e)* = A*^d^ + 3*^81 + 0^. 

3.1 Laplace transform 

Consider the Laplace transform of the linearized equations, 

(3.2) {L{£)-X)U = 0, AgC, xeR, U{e,x,X) e C^. 
Equation (3.2) can be cast as a first-order ordinary differential system in MJ, 

(3.3) W' = A{e,X)W, A G C, x G M, W{e,x,X) €C\ 
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where the Hmits A± of A at ±00 are given by 
(3.4) 












where |± denotes evaluation at U^, b± := b{U^), and J is defined in (1.4). 

Considered as an operator in L^(R;C^), L is closed, with domain iJ^ dense in 
L^. Similarly, for all A, the operator 



d_ 
dx 



- A(A) : H\R; C^) c L'^(R; C^) ^ L^{R; C^) 



is closed and densely defined. 

The following straightforward Lemma gives a correspondence between (3.2) and 
(3.3). 

Lemma 3.1. Let X e C and f = (/i, /2) G ^^(M; x C^). // the equation 
(3.5) (L -X)U = f 

has a solution U =: {r,w) G H^{R;C'^ x C^), then W := {T,w,bw') G H'^{R;C'^) 
satisfies 



(3.6) 



W' = A{X)W + F, 



with F = (/i,0,/2) G L2(M;C^). Conversely, let F = (/i,0,/2) G L'^(R;C'^) and 
X£C. IfW = {wi,W2) G H'^{R]C'^) satisfies (3.6), then a solution in if2(M;C^) to 
(3.5) mi/i / = (/i, 72) ^iwen 6j/ C/ = wi. 

Similarly, the dual eigenvalue equation 

(3.7) {L{e)*-X)U = 0, A G C, y G M, C/(£,y,A)GC^ 
can be cast as 

(3.8) W' = A{e,X)W, A G C, y G M, l^(£, y. A) G 
where the limits A± of A at ±00 are given by 



/ -s-'X 



(3.9) 







r>— 1 \ 7tr \ Q „tr xtr i „— 1 Ttra ^tr\i,tr— 1 



y s-'XJ^- X - d^g^l -{dyjf^± + s-'j^-drfll)bT' J 
A correspondence between (3.7) and (3.8) holds, as in Lemma 3.1. 
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3.1.1 The limiting, constant-coefficient equations 

Associated with (3.2) and (3.7) are the hmiting, constant-coefficient eigenvalue equa- 
tions 

(3.10) (L±(e) - X)U = 0, 
and 

(3.11) {L±{e)* - X)U = 0. 

Definition 3.2 (Normal modes). We call normal modes the solutions (A, C/) of 
equations (3.10) and dual normal modes the solutions (A, J7) of equations (3.11). 

Associated with (3.3) and (3.8) are the hmiting, constant-coefficient differential 
equations 

(3.12) W' = A±{£,X)W, 
and 

(3.13) W' = A±{e,\)W, 

where A± and A± are defined in (3.4) and (3.9). 

There is a correspondence between solutions of (3.10) and solutions of (3.12): 

Lemma 3.3. //(Ao,f/), U =: {t,w), is a normal mode, thenW := {T,w,bw') solves 
(3.12) at X = Aq. Conversely, if W = iwi,W2) G x solves (3.12) at X = Xq, 
then (Ao,'u^i) is a normal mode. In particular, 

(i) Eigenvalues ^ of A± satisfy 

(3.14) det(-/LtA± + i?B± + G±-X) = Q, 

and associated eigenvectors, satisfying A-|-(A)W^ = ^xW, have the form W = 
{U,W2) eC^x C3, with 

(3.15) U G ker(-(Li^± -|- /i^5± + G± - A), U =: (r, w), W2 := nb±w, 

(ii) Normal modes (A, U) satisfy 

(3.16) C/ = ^e^''^''(^)?7f (x,A), 

j 

where the are eigenvalues of A±, and the Up are polynomials in x. 



19 



The correspondence between (3.11) and (3.13) is similar. In particular, eigen- 
values jl of A± satisfy 

(3.17) det{fiA*^ + il^B*^ + - A) = 0, 

associated eigenvectors, satisfying K±{X)W = fiW , have the form W = {U,W2) G 
X C^, with 

(3.18) U G keiiflA*^ + ft^B^ + G*^- A), U =: (f, w), W2 := nb^^w, 
and dual normal modes satisfy 

(3.19) C7 = ^e^^^Wc7±(j/,A), 

j 

where the /i^ arc eigenvalues of A±, and the are polynomials in y. 

If jl{X) is an eigenvalue of A-i-(A), then /i(A) = —//(A), where /li(A) is some 
eigenvalue of A±(A). The matrices A±,B± and G± having real coefficients, the 
complex conjugate of is an eigenvalue of A-|-(A). We can thus relate the solutions 
of (3.14) and (3.17) by 

/i(A) = -^(A). 

Note that = 1, (t>(T^) = 0, and 0'(T^) = imply that the v derivative of the 
coupling reaction term k(l>(T)z vanishes when evaluated at U^. In particular, in 
{v,z) coordinates, 

with the notation of Section 1.3, \± denoting evaluation at Uj., and (j)± := 4>{Tj_), so 
that = 0, while by (1.8) and Assumption 1.2, > 0. This triangular structure 
of the matrix —iiA± + ij?B± + G± allows a simple description of the solutions of 
(3.14). Indeed, (3.14), a polynomial, degree four equation in A, splits into the linear 
equation 

(3.20) ixs + - fc^± - A = 0, 
and the degree three equation 

(3.21) del{-iidj\^ + i?h\^^ - A) = 0. 

By inspection, (3.20) is quadratic in while (3.21) is degree five in ji. Thus, the 
four solutions A(/x) of (3.14) correspond to seven eigenvalues /x(A) of A(A). 
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3.1.2 Low-frequency behaviour of the normal modes 

We describe here the behaviour of the normal modes in a smah bah B{0, r) := {A G 
C,|A|<r}. 

Definition 3.4 (Slow modes, fast modes). We call slow mode at ±00 any family of 
normal modes 

{{X,U{X)}x(zBio,r), for some r>0, 

such that, in (3.16), ^^(0) = 0, for all j. Normal modes which are not slow are 
called fast modes. We define similarly slow dual modes and fast dual modes, using 
(3.19). 

The solutions of (3.20) are 

(3.22) /x± = ^(-s + (s2 + 4d(A + fc0±))i), 

(3.23) 4 = -^is + (s' + 4d(A + A;0±))^); 

they depend analytically on A (in the case of and /j^, this is ensured by s > 0, 
assumed in Definition 1.3), and satisfy, for A in a neighborhood of the origin, 

(3.24) /i+ = s-^A - s-^dX^ + 0{X^), //J > 0, fif < 0. 

Note that the inequality /ij > is a consequence of ^_ > 0. By (3.18), the eigen- 
vector of A+ that is associated with — is 

(3.25) L+ = ( ) G X C', |+(0) = it, 
where 

(3.26) it ■■= {0 1 )*"" 

is the reactive left eigenvector of associated with the reactive eigenvalue of A^. 
We label LJ, Lf the eigenvectors of A± associated with —^^ and — /U^. By the 
block structure of — + ijl^B± + G±, spectral separation of /X4 and fi'^ (and of 
nt and //^), the eigenvectors and arc analytical in A, in a neighborhood of 
the origin (see for instance [Kat], II. 1.4); in particular, 

(3.27) tt = it + 0(A), ixth^ii = 0{\). 

The solutions of (3.21), seen as an equation in A, are the eigenvalues of the 
matrix —fidyf^_^ + ii^b\^^. By (1.3) and the block structure of A±, we find that the 

spectrum of dyf^^ is 

(^{dvfL) = {-s{£) - CT±, -s{e), -s{e) + a±}. 
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The eigenvalues of d^f^^ are distinct, hence, by Rouche's theorem, the eigenvalues 
of — + are analytic in /it, for small fi, with expansions 

Ai = s + cr± + f3ffi + 0(/^^), 
(3.28) X2 = s + pffx + 0{iJ^), 

A3 = s-a± + pfii + 0{n^). 

By (H3) (Section 1.5), > for all j. Inversion of these expansions yields analytic 
functions /i^, called fluid modes, and defined in a neighborhood of the origin in C;^ : 

nf := {s + a±y^X - (s + a±y^pt>''^ + 0{X^), 



± 



(3.29) /i^ 



By (3.18), the eigenvectors of A that are associated with these eigenvalues are 



= {s- a±r'X -{s- a^r^PiX" + 0{X^). 



(3.30) Lf{X) = ) G X C\ efiO) = if, 1 < i < 3, 

where the vectors if, if and if are the left eigenvectors of A± associated with the 
fluid eigenvalues —s — a±, —s, and —s + a±; they have the form 

(3.31) if :={ * * * )*^ 1 <i <3. 

The eigenvalues of —d^f^j^ + being distinct, the associated eigenvectors are 

analytical as well, so that the Lf, 1 < j < 3, are analytical in A; in particular, 

(3.32) if = if + 0{X), i^fb'^if = OiX). 

Finally, the equation det{— fid^, f^j^ + fi'^b\^_^) = has two non-zero solutions 
7^,7^, corresponding to the remaining two (fast) modes, solutions of 

(3.33) KT^'^c-^su/j? + (kc~^ (s^ - rT±^e±) + vs^^^^n + s{s'^ - cr^) = 0. 

The Lax condition (1.9) implies that solutions of (3.33) are distinct and have small 
frequency expansions 



^if = -if + 0(A), 7f > 0, 7+ < 0. 



Wc label Lf and hf the eigenvectors of A associated with —^xf and — /Xy . Again, 
by spectral separation, hf and Lf are analytical in A. 
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Lemma 3.5. For some r > 0, equations (3.13) have analytical bases of solutions in 
5(0, r), 



(3.35) ■■= {Vth<j<7, Vj^:=e-^''f^^^Lf{X), 

where the eigenvalues /i^ are given in (3.22), (3.23), (3.29), and (3.34) and the 
eigenvectors associated with the slow modes are given in (3.25), (3.27), (3.30) and 
(3.32). 

Proof. The above discussion describes analytical families /i^, L^, such that the 
vectors defined in (3.35) are analytical solutions of (3.13). For A 7^ 0, the 
eigenvalues /x^ are simple, so that the families define bases of equations (3.13). 
By inspection of the expansions at A = 0, the families 5^ define bases of equations 
(3.13) at A = as well. □ 

The above low-frequency expansions of the eigenvalues show that 

(i) Equation W' = A-(\)W has a 3-dimensional subspace of solutions associated 
with slow modes {fij , j = 1,2,3) and 4-dimensional subspace of solutions 
associated with fast modes {fx^,fx^, iiQ,ix^); 

(ii) Equation W' = A+(A)W^ has a 4-dimensional subspace of solutions associated 
with slow modes (/x^, j = 1,2,3, and //J) and a 3-dimensional subspace of 
solutions associated with fast modes (/i^,/i^, /Xy"). 

3.1.3 Description of the essential spectrum 

We adopt Henry's definition of the essential spectrum [He]: 

Definition 3.6 (Essential spectrum). Let B he a Banach space and T : D{T) C 
B B a closed, densely defined operator. The essential spectrum of T, denoted by 
a ess (T) , is defined as the complement of the set of all A such that A is either in the 
resolvent set of T, or is an an eigenvalue with finite multiplicity that is isolated in 
the spectrum ofT. 

By Lemma 3.3, the matrix A-i-(A) has a non trivial center subspace if and only 
if AgC±, 

C± := {A G C, det{-i^A± - fB± + G±-\) = Q, for some ^ G R}. 

The following Lemma can be found in [He] (Theorem A. 2, Chapter 5 of [He], 
based on Theorem 5.1, Chapter 1 of [GK]): 

Lemma 3.7. The connected component of C \ (C- U C+) containing real +00 is a 
connected component of the complement of the essential spectrum of L{e). 
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Figure 3: Spectrum of L{e). 

The reactive eigenvalues of —i^A± — ^^B± + G± are 

X = i^s - $,^d - k(f)±. 
For small |^|, the fluid eigenvalues satisfy 

X = ia^- PC'^ + 0{f), aeR, P > 0, 
as described in Section 3.1.2; for large |^|, they satisfy 

(3.36) A = — ^^(a + 0{^~^)) (parabolic eigenvalues), 

with a G {vT^^, kc~^t^^}, or 

(3.37) A = is^ + 0(1) (hyperbolic eigenvalue). 

This implies that the essential spectrum is confined to the shaded area in Figure 3, 
the boundary of which is the union of an arc of parabola and two half-lines. (The 
origin A = is an eigenvalue, associated with eigenfunction (U^)'; the existence of 
bifurcation eigenvalues 7(e) ±iT(£) is assumed in Theorem 1.14, the proof of which 
is given in Section 5.) 

Remark 3.8. The essential spectrum, as given by Definition 3.6, is not stable un- 
der relatively compact perturbations (see [EE], Chapter 4, Example 2.2); namely, a 
domain of the complement of the essential spectrum of a ( closed, densely defined) 
operator T is either a subset of the complement of the essential spectrum ofT + S, 
or is filled with point spectrum ofT + S, where S is a relatively compact perturbation 
ofT. 
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Remark 3.9. By the Frechet-Kolmogorov theorem, L is a relatively compact pertur- 
bation of L± . ( This observaMon is the first step of the proof of Lemma 3.7, see Henry 
[He] J The pathology described in Remark 3.8 does not occur in the right half-plane 
here, as we know by an energy estimate that if X is large and real, A ^ ap{L). 

3.1.4 Gap Lemma and dual basis 

Let A be the connected component of C \ (C_ U C+) containing real +oo. 

Definition 3.10 (Stable and unstable subspaces at ±oo). Given A G AuB(0, r), r as 

in Lemma 3.5, denote by S{A±{X)) the stable subspace ofA±{X) (i.e., the subspace of 
generalized eigenvectors associated with eigenvalues with negative real parts ) and by 
U{A±{X)) the unstable subspace of A± (i.e., the subspace of generalized eigenvectors 
associated with eigenvalues with positive real parts). We define similarly 5'(A-|-(A)) 
and C/(A±(A)). 

By definition of C±, given A G A, the matrices A-|-(A) do not have purely imagi- 
nary eigenvalues, so that S'(A±(A)) ® ?7(A±(A)) = C^, and S'(A±(A)) ® C/(A±(A)) = 
for ah A G A. 

Lemma 3.11. The vector spaces S{A±{X)) and U{A±{X)) have analytical bases in 
A. 

Proof. By simple-connectedness of A, the Lemma follows from a result of Kato 
([Kat], IL4), that uses spectral separation in A. □ 

Corollary 3.12. Equations (3.13) have analytical bases of solutions in A. 

Proof. Basis elements of the stable and unstable spaces defined in Definition 3.10 
are associated, through the flow of (3.13), with bases of solutions of (3.13). The 
matrices A± depending analytically on A, the flow of (3.13) is analytical in A. □ 

Lemma 3.13. For X real and large, dimS'(A_|_(A)) = dimS'(A_(A)) = 3. 

Proof. Prom Lemma 3.3, fi is an eigenvalue of A±(A) if and only if A is an eigenvalue 
of —ijlA± -\- n^B± -\- G±. As in Section 3.1.3, for large the eigenvalues of —fiA± -\- 
^PB± + G± are s/i + 0(1) (hyperbolic mode) and i/t^^/j? + 0{^), KT^^c~^fj,'^ + 
0{jj,), dfj? + 0{ii) (parabolic modes), c~^ as in Assumption 1.1. Inversion of these 
expansions gives three stable eigenvalues for both A_ and A+. □ 

Remark 3.14. The above Lemma implies in particular that A is a domain of con- 
sistent splitting, as defined in [AG J]. (See also Section 3.1 of [LyZ2].j 

Given A G A, the flow of (3.13) associates basis elements of S'(A_)_(A)) with solu- 
tions of (3.13) which are exponentially decaying as t — > -|-oo, and basis elements of 
?7(A_(A) with solutions which are exponentially decaying as t ^ — oo. Similarly, the 
spaces 5(A_(A)) and t/(A+(A) are associated with exponentially growing solutions, 
at — oo and +oo respectively. 
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Figure 4: Normal modes on the — oo side. 



Definition 3.15 (Decaying and growing normal modes). We call decaying dual 
normal mode at ±00 any continuous family of dual normal modes {A, U{X)}, A € 
B{0,r), r as in Lemma 3.5, such that for a// A G A n B(0,r), U(X) corresponds 
to a decaying solution of (3.13) at ±00. Families of normal modes which are not 
decaying are growing. We define similarly decaying dual normal modes and growing 
dual normal modes. 

By continuity of the eigenvalues and spectral separation in A, if for some A G A 
a continuous family of normal modes corresponds to a decaying (resp. growing) 
solution, then it corresponds for all A G A to a decaying (resp. growing) solution. 

By (1.9), (3.24) and (3.29), fi^, /Xg", and /xj are growing (in the sense of 
Definition 3.15) at +00, while /Xg" and fi^ are decaying. 

Similarly, ^J, and ^Iq are growing, while , A*4 ^^'^ 1^7 ^^e decaying. 

The normal modes with which the characteristics of (1.1) are associated are 
pictured on Figures 5 and 4. In particular, slow normal modes associated with 
incoming characteristics are growing. 

Definition 3.16 (Normal residuals). A map (y, A) — >■ &'^{y,X) G defined on 
[yo, +00) X B{0, r), for some yo > 0, r > 0, is said to belong to the class of normal 
residuals if it satisfies the estimates 

|e+|<C, l^j^G+l <C(|A|+e-^l2'l). 

for some 9 > and C > 0, uniformly in y > yo and A G B{0,r). 

We define similarly the class of normal residuals on (—00, —yo) x B[Q,r). 
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Figure 5: Normal modes on the +00 side. 



Lemma 3.17 (Fast dual modes). Equations (3.8) has solutions 

, , Wq , ( growing ) and , Wq , Wy" ( decaying ), 
which for A G B(0,r), r possibly smaller than in Lemma 3.5, satisfy 

(3.38) #f = e-y^f^'^LfiO) + e-^l^le±. + AGJ), V ^ ±yo, 

for some yo > independent of X, where the constant vectors L^i^) are defined in 
Section 3.1.2, and 0^, @fj are normal residuals in the sense of Definition 3.16. 

Proof. With the description of the normal modes in Lemma 3.5, this is a direct appli- 
cation of the Gap Lemma (for instance in the form of Proposition 9.1 of [MaZ3]). □ 

Lemma 3.18 (Slow dual modes). Equation (3.8) has solutions 

Wi , W2 (growing) and W3" , , W2 , , ( decaying ) , 
which for X G B{0,r), r possibly smaller than in Lemma 3.5, satisfy 

(3.39) T4f = e-^f^f^^^ (L±(0) + Xef) , y ^ ±yo, 

for some yo > independent of X, where the constant vectors L^(0) are defined in 
Section 3.1.2, and Qf are normal residuals. 
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Proof. The Conjugation Lemma ([MeZ]; Lemma 3.1 of [MaZ3]) implies that there 

exists a family of matrix-valued applications {®~^ ^)}\^B{o,r)7 foi" some r > 
possibly smaller than in Lemma 3.5, such that the matrix Id + 0+ is invertible for 
all A and y, the application 0+ is smooth in y and analytical in A, with exponential 
bounds 

\d{d^@^\ < CjkC-'^y, for some 6 > 0, Cjk > 0, for y > yo, 
for some yo > 0, and such that any solution W of (3.8) has the form 

(3.40) = (Id + @+)V+, for y > yo, 

where is a dual normal mode, and, conversely, if is a dual normal mode, 
then W defined by (3.40) solves (3.8) on y > yo. 

Equation (3.8) at A = has a four-dimensional subspace of constant solutions; 
let {Wj}i<j<4 be a generating family. The normal modes with which, through 

(3.40) , the Wj are associated are slow normal modes. Hence, by Lemma 3.5, there 
exist coordinates Cjk such that 

W° = (Id + 0+(-,O)) cjkL+iO), y>yo, 

l<k<4 

which implies in particular that the matrix c := (cjfc)i<i,fc<4 is invertible. Then, for 

1 < i < 4, 

(Id + ©+(-,0))L;(0)= {c-')jkW^, 

l<fc<4 

in particular, (Id + ©"'"(•, 0))L^(0) is constant, hence, by exponential decay of 
equal to L'j'{0). We can conclude that, for 1 < j < 4, 

w+ := (Id + ©+)v;+ 

(where V^^ is defined in Lemma 3.5) is a solution of (3.8) on y > yo; which can be 
put in the form (3.39). 

The proof on the — oo side is based similarly on the decomposition of the fluid 
components of the Wj onto the (fluid) dual slow modes Vj~ , for 1 < j < 3. □ 

3.1.5 Duality relation and forward basis 

We use the duality relation, introduced in [MaZ3] , 

(3.41) W^^'SW = 1 
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that relates solutions W of the forward equation (3.3) with solutions W of the 
adjoint equation (3.8) through the conjugation matrix in (r, w, bw') coordinates 





' -An 


-A12 







-A21 


-A22 


Idc3 




V 


-Idc3 






where A is the convection matrix defined in (1.16). Namely, is a solution of (3.3) 
if and only if it satisfies (3.41) for all solution W of (3.8), and conversely is a 
solution of (3.8) if and only if it satisfies (3.41) for all solution W of (3.3). (See 
Lemma 4.2, [MaZ3]; note that the reactive term contains no derivative, hence does 
not play any role here.) 

Remark that there exist vectors such that 

(3.42) ^tA±r^ = -Sjk, l<j,k<A. 
Let be vectors of the form 

(3.43) := e-'@f„ 

where for 1 < A; < 4, are given by (3.43), and where 9^^. are normal residuals. 
With the notation of Lemmas 3.17 and 3.18, let 



(3.44) L 



^ ( Lf{0) if is slow. 



[ Lf{0) + e-^l^'ej. if fif is fast. 
Lemma 3.19 (Forward and dual basis). For some r > and yo > 0, 

• equation (3.3) has analytical bases of solutions {W^, . . . ,W'^}x^AuB{o,r)y for 

y ^ ±yo; 

• equation (3.8) has analytical bases of solutions {W^ , . . . ,W^}xeMJB(o,r): for 

y ^ ±yo, 

such that for A G B{Q,r), 

(3.45) =e^^f^^\Rf + Xef), y^±yo, 

(3.46) =e-y^f^^\Lf + Xef), y^±yo, 



where and are defined in (3.43) and (3.44), and @^ and are normal 



residuals; the fast forward modes and satisfy also 

(3.47) W^ix, A) = ( ^^'l^'^^ ) + XQfix, A), x ^ ±yo, 

where \ef\ + \d^ef\ < Ce"'^'^!, for some C,9 > 0, uniformly in X € B{Q, r). 
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Proof. Given a family {Fx, . . . , F-j} of vectors in C^, let co\{Fj) denote the 7x7 
matrix col(Fj) := ( Fi ... Fj ) . 

Let and Wp as in Lemma 3.17 and 3.18. For all A € AUi?(0, r), the families 
. . . , W~} and {M^j^, . . . , W^n*"} are bases of solutions of (3.8), on y < — yo and 
y ^yo respectively. In particular, the 7x7 matrices W'^='= := col(W"j^) are invertible 
for all A G -6(0, r) and y ^ ±yo. Let 

(3.48) := ((WO±)*^5)-i =: co\{Wj^^). 

For the forward modes Wj^ defined in (3.48) to satisfy the low-frequency description 

(3.49) Wf^ = e^^f^^\R^^ + e-'\-\@lf + X@lf), y ^ ±yo, 

where -R°^ are constant vectors and 6°^ are normal residuals, it suffices, by (3.41), 
that the matrices R°=^ := col(-R°='=) and ©^"^ := col(e°=^) satisfy 

(3.50) L*^5R° = Idc7, 

(3.51) (L + e-^l^lei)*^50? = -0f5R°, 

(3.52) (L + e-'^^^le? + Ae2)*"50^ = -ef5R°, 

where L± := col(L=^(0)) and G)± := col(e^^.) appear in the low- frequency description 

of the Wp. In (3.50)-(3.52), the ± exponents are omitted. The matrices being 
invertible, (3.50) (with -|- or — ) has a unique solution, and, for yo large enough and 
r small enough, equations (3.51) and (3.52) have unique solutions in the class of 
normal residuals. Note that for 1 < j,k < 4, equation (3.50) reduces to (3.42), up 
to exponentially decaying terms, so that the vectors have the form (3.43). 

Remark now that {U^Y satisfies L{e){U^)' = 0, and decays at both — oo and 
-l-co, hence {U^)' is associated with decaying fast normal modes; by Lemma 3.17, 
there exist constants c^, such that 

(3.53) ('^T'")=''"'^°"*-= ^'^^'^^'l-"- 

We may assume, without loss of generality, that ^ 0. Let now 

W" := ( W^- W^- clWl' W^- W^- W^- ), 

W+ := ( <+ Wl+ <+ <+ <+ El=5c;W^°+ ), 

and W='= =: co\{Wp). These forward modes satisfy (3.45) and (3.47). Let finally 
W^*^ := (5W^)^^ =: co\{Wp), so that, in particular, the slow modes of and 
W"^ coincide. Wc can prove as above that the low-frequency description (3.49) of 
the forward modes carries over to the dual modes through the duality relation, so 
that (3.46) is satisfied. □ 
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3.1.6 The resolvent kernel 

Let 

L'^{n, V'{R)) := {(p £ V\n X M), for all cp G V{R), (0, ^p) G L'^{n)} 

A linear continuous operator T : L'^{R) ^ ^^(R) operates on ^^(M, by 

r(,)( 



(T^,V?) := r(^,(/?). Let Tr.)6 G L2(M,X)'(M)) be defined by {Tj,ip) = ip{x), for all 



X G JR. 

Definition 3.20 (Resolvent kernel). Given X in the resolvent set of L{e), define the 
resolvent kernel Q\ of L{e) as an element of L'^{Ra:,'D'{Ry)) by 

g^:= iL{s)-Xr\)S. 

Given y G M, let Sy = sgn(y), and 

D(y) := {j, II^ slow and decaying}, 

so that 

D(y) = {3}, ify<0, D(y) = {1, 2, 3, 4}, if y > 0. 

Given x, y G M, let D(x, y) be the set of all (j. A;) such that for all x, y, for > 
and |A| small enough, SR(/Lt|^x — /x^^y) < 0, that is, 

D(x,y) = {(j, fc), /i^^ and slow and decaying} 

IJ {(j,i), = Sy, |y| < l-^l, /"f slow and decaying} 
lj{(j>j), Sa; = |x| < |y|, /x^^ slow and decaying}, 

so that 

{ {(1,1), (2,2), (1,3), (2,3)}, x<y<0, 
{(1,3), (2, 3), (3, 3)}, y<x<0, 
0, y<0<x, 
{(j,A:), l<j<4,l<A;<2}, x < < y, 
{(j,j), l<j<4}, 0<x<y, 
V 0, 0<y<x. 
Define now the excited term 

5,(x,y):=A-HC/^)'(x) ^ [c^.J^^^^^e"^'^? W, 

ieD(2/) 

and the scattered term 

(i,fc)GD(x,y) 

where the vectors £^ arc defined in (3.26) and (3.31), the vectors r^ arc defined in 
(3.42), and the transmission coefficients [c^ _|_] and [c4'^] are constants. 



D(x,y) := < 
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Proposition 3.21. Under (1.18), /or A G B{0,r), the radius r being possibly smaller 
than in Lemma 3.5, there exist transmission coefficients [c^ _|_] and [<4'^] such that 
the resolvent kernel decomposes as 

where TZ\ satisfies 

\dxdynx\ < Ce-^l^-^'l + CA"'e-^l^l ^ e"*"^? 

{j,k)er){x,y) 

for a G {0, 1, 2, }, a' G {0, 1}, for some C,0 > 0, uniformly in x, y and A G -B(0, r). 

Proof. The duality relation (3.41) allows to apply Proposition 4.6 of [MaZ3] (and its 
Corollary 4.7), which describes Gx as sums of pairings of forward and dual modes, 
for A in the intersection of A and the resolvent set of L. By Lemma 3.19, Qx extends 
as a meromorphic map on 5(0, r). 

The excited term £x comprises the pole terms, corresponding to pairings of a 
fast, decaying forward mode associated with the derivative of the background wave 
with a slow, decaying dual mode, i.e. /W^ for y < and /W^' for y > 0, 
1 < j < 4. 

The next-to-leading order term is the scattered term S\. It corresponds to pair- 
ings of a slow forward mode with a slow dual mode. For y < 0, the scattered term 

comprises only fluid modes. For y < < .t and for < y < .x, the scattered term 
vanishes, as there are no outgoing modes to the right of the shock (see Figures 1 
and 4). 

By the Evans function condition (1.18) and Lemma 6.11 of [MaZ3], the residual 
TZ\ does not contain any pole term; it comprises: 

(a) the contribution of the normal residuals to the fast forward/slow dual pairings 
involving the derivative of the background profile, 

(b) the fast forward/slow dual pairings not involving the derivative of the back- 
ground profile, 

(c) the contribution of the normal residuals to the slow forward/slow dual pairings, 
and 

(d) the slow forward/fast dual pairings. 

Term (a) is bounded by the first two terms in the upper bound for TZx. Term (b) 
is smaller than term (a) by a 0(A) factor. Term (c) is bounded by the third term in 
the upper bound. By the Lax condition (1.9), the Evans function condition (1.18) 
and Lemma 6.11 of [MaZ3], term (d) is also bounded by the third term. □ 
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3.1.7 The Evans function 

By Lemma 3.13, for all A G A, the dimensions of i7(A_(A)) and S'(A+(A)), the vector 
spaces associated with decaying solutions of (3.3) at — oo and +oo, add up to the 
full dimension of the ambient space: 

dimC/(A_(A)) + dim5(A+(A)) = 7. 

Definition 3.22 (Evans function). On A U B{0,r), define the Evans function as 

Die, A) := det(W^f , W^", W4-, T^f, ^^5+, W^, W^+)^,=o. 

The Evans function D satisfies Proposition 1.8; it has a zero at A = 0, as reflected 
in equality (3.53). 

3.2 Inverse Laplace transform 

Similarly as in Section 3.1.6 (or Section 2 of [MaZ3]), define the Green function of 

L{£) as 

(3.54) g := e*^(")r(.)<5, 

where {e*^^^^}t>o is the semi-group generated by L{£). That is, the kernel of the 
integral operator 6*"^^^^ is the Green function Q evaluated at t = tg- 

Assuming (1.18), the inverse Laplace transform representation of the semi-group 
by the resolvent operator (see for instance [Pa] Theorem 7.7; [Z3] Proposition 6.24) 
yields 

(3.55) g{e, X, t- y) = ^P-V. / e^'Gxie, x, y) dX, 
for 770 > sufficiently large. 

3.2.1 Pointwise Green function bounds 

Introduce the notations 

[y _ 2 

errfn(y) := / e ^ dz, 
J —00 

and let, for y < : 




i<i<4 V \\/WiJ Vv^/^i* 
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and 

(3.58) £{e,x,t;y) := {U')' {x)e{e,t;y). 

In (3.56)-(3.57) and below, the {o-^}i<j<4 are the characteristic speeds, i.e. the hmits 
at ±00 of the eigenvalues of djjF{e, Lf^), ordered as in (1.3), the P^, 1 < j < 3, are 
the positive diffusion rates that were introduced in (3.28) (and which depend on e, 
as do the characteristic speeds), and (3^ := d, the species diffusion coefficient. 

Let for y < : 
(3.59) 

«5 := X{t>i} E r7£T-(4vr/37t)-|e-(--^-'^7*)V4/37t 

l<j<2 

+ E [C]-7^r^(4^/?^;:o-^e-(--^;=)V4/3^;i*), 

i<i<2 

and for y > : 
(3.60) 

^■■=Xi.>n^^ E r/f;"(4<«)-ie-(—r<)'/4,i;. 
l<j<4 

+ ^^te'ipS? E E [4';i.-.-<r'(''-'5j:';r4e-('-*;iv<?., 

l<j<4 l<k<2 

where the indicator function X{t>i} is identically equal to 1 for t > 1 and otherwise, 
and 

Let 

(3.61) H := h{s, t, X, y)T^+stS, hit = 0> 

where the notation r(.)(5 was introduced at the beginning of Section 3.1.6. 

Let finally So be the scattered term defined in (3.59)- (3.60) in which [c^'^] = 1 
for all j, k. 

Proposition 3.23. Under (1.18), there exists transmission coefficients [c^ ±] and 
[4',±]' satisfying 

(3.62) J rl = [cO,J« - ^;i) + \c\-\r^ + 

. = [cj+l = 0, 
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where l<k<3if-k = + and l<k<2if-k=—, Uj. =: (f^;2;j_), such that 
the Green function Q{e,x,t;y) defined in (3.54) may be decomposed as a sum of 
hyperbolic, excited, scattered, and residual terms, as follows: 

(3.63) g = n + s + s + n, 

where Ti., 6 and S are defined in (3.56)-(3.61), with the estimates 

\d^d^dfh\<Ce-'\ 

+ C ((r^ + e-^l-l)i-^(°+"') + a't-h-^\y\) \So\. 

uniformly in e, for k G {0, 1}, a G {0, 1, 2}, a' G {0, 1}, for some 6,C,M > 0. 

Proof. We only check (3.62), as decomposition (3.63) and bounds (3.64) are easily 
deduced from Proposition 7.1 of [MaZ3] and Proposition 7.3 of [LRTZ]. (See also 
Proposition 3.7 of [TZ2], especially equations (3.30)-(3.33) and (3.38).) 

The description of the residue of Qx aX \ = fov y < and y > implies 

1<J<3 

corresponding to equation (1.34) in [MaZ3]. The (reactive) left eigenvector vector 
£4" being orthogonal to the (fluid) left eigenspace span{^^}i<j<3 (see (3.26) and 
(3.31)), this implies (3.62) (i). 

Given Uq & L^, the estimates for H and TZ imply 

lim / {n + n)Uodydx = 0. 

Hence, by conservation of mass in the fluid variables, (1.5), for all Uo £ L^, 

(3.65) lim / Tr{£ + S)Uo dy dx = / irUody. 

where tt : ^ ^ is defined by 7r{v,z) := v. (Equation (3.65) corresponds to 
(1.33) and (7.60) in [MaZ3].) Taking Uo G span{£=^}i<j<3, we find (3.62)(ii), and 
taking Uq parallel to i^, we find (3.62) (iii). □ 

Remark 3.24. The terms £ and S correspond to the low-frequency part of repre- 
sentation of Q by inverse Laplace transform of the resolvent kernel Q\, while the 
term Ti corresponds to the high-frequency part. As observed in [MaZ3], for low fre- 
quencies, the resolvent kernel in the case of real (physical) viscosity obeys essentially 
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the same description as in the artificial (Laplacian) viscosity case, hence the esti- 
mates on £ and S follow by the analysis in [LRTZ] of the corresponding artificial 
viscosity system, specialized to the case of strong detonations ( more general waves 
were treated in [LRTZ] The estimate of the terms H and TZ follows exactly as for 
the nonreactive case treated in [MaZ3, Z2]. 

Remark 3.25. Bound (3.64)(ii) is implied by bounds (7.1)-(7.4) of Proposition 7.1 

o/[MaZ3] and bounds (3.30), (3.32) and (3.38) of Proposition 3.7 of [TZ2]. 

Here the contribution of the hyperbolic, delta-function terms to the upper bounds 
for the spatial derivatives of TZ is absorbed in H, and the short-time, t < \a^\\y\, 
contributions of the scattered terms are absorbed in the generic parabolic residual 
term e-^te-^^-v? m _ 

Corollary 3.26. The excited terms £\ and £ contain only fluid terms: £^£2 = 
and££2 = 0. 

Proof. The equality ££'1 = follows from (3.62) (i). The resolvent kernel Qx is 
the Laplace transform of the Green function Q, so that the coefficients [c^ ±\ in 
Propositions 3.21 and 3.23 must agree. Hence. (3.62) (i) implies also £x£4 =0. □ 

Corollary 3.27. For all 7] > 0, for some C,M > 0, some 9i{r), s) > 0, the following 
bounds hold, for a £ {0, 1, 2} : 

Proof. By (3.62) (iii), the contribution of the reactive modes to S is 
(3.67) ;^^^^,^;^^^^,^_^r+4-(47rdt)-V2e-(-^+^*)V4<it. 

Given < x < y, we can bound e~'?2'e""(^~2/+st)^/4<it^ foj. \x — y\ < ^st, by 
and, for \x — y\ > ^st, by 

g-TO < g-»?l2/-^l/2g-r)j//2 ^ g-?7|j/-a;|/2g-?7st/4^ 

and this implies (3.66) (i). To prove (3.66) (ii), we note that the contribution of 
the parabolic terms t-2{t + l)-2e-(^-2')'/^* and Sq to 7^£| comes from Riemann 
saddle-point estimates of the sole scattered terms Sxl^ (see the proof of Proposition 
7.1 in [MaZ3] for more details). Hence (3.66) (i) implies (3.66) (ii). □ 

Remark 3.28. The proof of Proposition 7.1 o/ [MaZ3] shows that Proposition 3.23 
applies more generally to linear operators of the form (1.15) that satisfy (1.18) and 
the conditions (Al)-(A2), (H1)-(H4) of Section 1.5. 
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3.2.2 Convolution bounds 



Prom the pointwise bounds of Proposition 3.23 and Remarks 3.26 and 3.27, we 
obtain by standard convolution bounds the fohowing — > estimates, exactly as 
described in [MaZl, MaZ2, MaZ3, MaZ4, Z2] for the viscous shock case. 

Corollary 3.29. Under (1.18), for all t > 1, some C > 0, any rj > 0, for any 
l<q<P,l<P< +00, and / G L"? n W^'P, 



I 

JM. 



{S + n){;t;y)f{y)dy 

JM. 

[ dy{S + n)i;t-y)f{y)dy 
iS + n){;t;y)eifiy)e-'yUy 
ni;t;y)fiy)dy 



/ 

JR 



LP 



LP 



LP 



LP 



< Ct--^^-^--p^~~^\f\L. + Ce-^V\Lp, 



1/1 1 \ 1 



< Ct-^^-^-p-^\f\L. + Ce-'^V\Lp, 

< ce-^y\LP, 



where y'^ := max(y, 0). Likewise, for all x and all t >0, 

\dye{;t)\LP + \dte{;t)\LP < Ct-^^^^-p\ 

\dtdye{;t)\LP < Cr^^'-p^-^. 



4 Stability: Proof of Theorem 1.12 

We often omit to indicate dependence on e in the proof below. All the estimates 
are uniform in e. 



4.1 Linearized stability criterion 

Proof of Theorem 1.12: linear case. Sufficiency of (1.18) for linearized orbital sta- 
bility follows immediately by the bounds of Corollary 3.29, exactly as in the viscous 
shock case, setting 

6{t) := / e{x,t;y)Uo{y)dy 

JR 

SO that 

U-5{t)U'= [ {n + S + n){x,t-y)Uo{y)dy- 

JR 

see [ZH, MaZS, Z2] for further details. Necessity follows from more general spectral 
considerations not requiring the detailed bounds of Proposition 3.23; see the discus- 
sion of effective spectrum in [ZH, MaZ3, Z2]. The argument goes again exactly as 
in the viscous shock case. □ 
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4.2 Auxiliciry energy estimate 

Consider U the solution of (1.2) issued from Uq, and let 

(4.1) U{x,t) ■.= U{x + 5{t),t) -U{x). 
Then, the following auxiliary energy estimate holds. 

Lemma 4.1 (Proposition 4.15, [Z2]). Under the hypotheses of Theorem 1.12, as- 
sume that Uq £ H^, and suppose that, for < t < T, the suprema of \5\ and the 
norm of U each remain bounded by a sufficiently small constant. Then, for all 
0<t<T, for some 9 > 0, 

\U{t)\l, < Ce-''\Um'HS+C f e-'^'-'\\U\l. + \5\'){s)ds. 

Jo 

4.3 Nonlinecir stability 

Proof of Theorem 1.12: nonlinear case. Let U be the perturbation variable associ- 
ated with solution U as in (4.1); by a Taylor expansion, U solves the perturbation 
equation 

dtU -LU = d,Qf{U, d,U) + QriU) + 6{t){U' + d^cU), 
where the linear operator L is defined in (1.15), and 

(4.2) \Qf\<C\U\{\U\ + \d,U\), 
where C depends on ||C/||loo and ||t7||^yi,oo. 

Lemma 4.2. Under Assumptions 1.1, 1.2, if the temperature T associated with 
solution U satisfies \\T\\i^oo < Tj — r+ (by Lemma 1.4, < Ti — T^), then the 
nonlinear reactive term Qr has the form 

(4.3) Qr{U)=iie~'^"^^qr{U), 

where x+ := max(a;,0), ijq > is as in Corollary 1.6, and qr{U) = qr{w,z) is a 
scalar such that 

(4.4) \qr{U)\ < C\U^2, 
where C depends on o,nd \\U\\l°°. 

Proof. We use the specific form —k(t){T)z£'l of the reactive source in (1.1), together 
with Taylor expansion 

(0(r + T){z + z)- {(t>{f)z-{4>\f)T-z + 4>{f)z) 

= (p'{f)Tz + (p"{f + f3T)T'^z, 

for some < /3 < 1, and the fact that <t>'{f + T) < Ce-''o^"^ for |r| < Tj - T+, for 
> as in Corollary 1.6, by <p{T^) = together with the property that (p'iT) = 
for T <Ti and exponential convergence of U (x) to U^ as x ^ +oo. □ 
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Recalling the standard fact that U' is a stationary solution of the linearized 
equations (3.1), LU' = 0, or 

/ g{x,t;y)U'iy)dy = e*^C7'(x) = U'{x), 
Jr 

we have by Duhamel's principle: 

U{x,t)=5{t)U'{x)+ [ g{x,t;y)Uo{y)dy 

t 

+ [ [ g{x,t-s;y)i+e-^y^qr{U){y,s)dyds 
Jo, Jr 



* / dyg{x,t-s;y){Qf{U,d:,U) + 6U)iy,s)dyds. 
Jr 



Defining 

5{t) = - [ e{y,t)Uo{y)dy 
Jr 

(4.5) -/ [ eiy,t-s)ite-'^y^qr{U){y,s)dyds 

Jo^ Jr 

+ f [ dyeiy,t - s){Qf{U,d^U) + 6U){y,s)dyds, 
Jo Jr 

following [Z3, MaZl, MaZ2, MaZ4], and recalling Proposition 3.23, we obtain finally 
the reduced equations: 

Uix,t)= f {H + S + n){x,t-y)Uo{y)dy 

t 

' n{x,t- s;y){dy{Qf{U,d^U) + SU)+e+e-'^y^qr{U)) dyds 

I [{S + n){x,t- s; y)ete~''y^qr{U) dy ds 
Jo^ Jr 

- f [ dy{S + n){x,t-s;y){Qf{U,d^U) + 6U)dyds, 
Jo Jr 

and, differentiating (4.5) with respect to t, and recalling Corollary 3.26: 

m = - I dte{y,t)Uo{y)dy 
Jr 

+ [ [ dydte{y,t-s){Qf{U,d^U) + 5U){y,s)dyds. 
Jo Jr 



(4.6) 



(4.7) 



(4.8) 



5{t) = ~ [ e{y,t)Uo{y)dy 
Jr 

+ 11 dye{y,t-s){Qf{U,d^U) + 6U){y,s)dyds, 
Jo Jr 
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Define 

at) := sup ( \U{; s)\L.il + s)'^^'-'^^ + \6{s)\{l + s)U \6is)\) . 

0<s<t,2<p<oo ^ ^ 

We shall establish: 

Claim. There exists cq > 0, such that, for all t > for which a solution exists 
with uniformly bounded by some fixed, sufficiently small constant, there holds 

CW<co(|C/o|Lini/3+CW')- 
Prom this result, it follows by continuous induction that, provided 

\Uo\LinH^ < ^Cq, 

there holds 

(4.9) Cit)<2co\Uo\L^nH^ 

for all t > such that C remains small. For, by standard short-time theory/local 
well-posedness in iJ^, and the standard principle of continuation, there exists a 
solution U G on the open time-interval for which \U\fj3 remains bounded, and 
on this interval is well-defined and continuous. Now, let [0,T) be the maximal 
interval on which |f/|/f3 remains strictly bounded by some fixed, sufficiently small 
constant 6 > 0. By Lemma 4.1, we have 

mt)\ls < C\Umlse-'' + C f e-'^'-'H\U\h + \Sf){r)dT 

Jo 

<c'{\umH^+atf)ii+t)-K 

for some C,C',0 > 0, and so the solution continues so long as ( remains small, 
with bound (4.9), at once yielding existence and the claimed sharp U^r\H^ bounds, 

2 < p < GO. 

Proof of Claim. We must show that each of the quantities 
|f/|Lp(l + s)^^^"^\ |<j|(l + s)5, and \6\ 

is separately bounded by C(|[/o|Lin/i'3 + C(0^)) some C > 0, all < ,s < t, so 
long as C, remains sufficiently small. By (4.6)-(4.7) and the triangle inequality, we 
have 

\U\lv < la + lb + Ic + Id, 

< Ila + IIb, 
\S{t)\ < Illa + IIIb, 
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where la is the LP norm of the first integral term in the right-hand side of (4.6), 
lb the second term, etc., and similarly Ila is the modulus of the first term in the 
right-hand side of (4.7), etc. 

We estimate each term in turn, following the approach of [MaZl, MaZ4]. 

The linear term satisfies bound 

Ia<C|^o|LinLp(l + i)"^^'"^\ 

by Proposition 3.23 and Corollary 3.29. 

Likewise, applying the bounds of Corollary 3.29, we have 

lb < CC{tf f e-''(*-^) (1 + s)-^ds 
Jo 

<catf{i+trK 

and (taking g = 2 in the second estimate of Corollary 3.29) 

Ic + Id < C f e-i^^-'\\U\Loo + \d^U\Loo + \5\)\U\Lv{s)ds 
Jo 

+ C [ {t- s)-t+i(|C/|^oc + \5\)\U\Hiis)ds 
Jo 

<CC{tf [ e-''(*-')(l + s)"^(^"p)~^ds 
Jo 

+ CC{tf [ (t-s)"^+^(l + s)-3ds 
Jo 

<ccit)\i + t)-'^^^-p\ 

na< \dte{t)\Loo\Uo\L^ < C\Uo\L^{l+t)--2 

and ^ 

nb< / \dydte{t-s)\L2{\U\Lo<^ + \5\)\U\Hi{s)ds 
Jo 

< CCitf [ (t-s)-t(l + s)-!(is 
Jo 

<cc{tf{i+t)-K 

while 

nia< \e{t)\L^\Uo\Li < C\Uo\l1 

and ^ 

IIIb< / \dye{t-s)\L2{\U\L^ + \6\)\U\Hi{s)ds 
Jo 

ft 1 3 

<CC{tf {t- s)~{l + s)-^ds 
Jo 
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This completes the proof of the claim, establishing (1-19) for p > 2. The remain- 
ing bounds 1 < p < 2 then follow by a bootstrap argument as described in [Z2]; we 
omit the details. □ 



5 Bifurcation: Proof of Theorem 1.14 

Given two Banach spaces X and Y, we denote by C{X, Y) the space of linear contin- 
uous applications from X to Y, and let C{X) := C{X,X). We use (1.20) to denote 
weighted Sobolev spaces and norms. Let a;"*" := max(0,x). Given a constant 77 > 
and a weight function a; > 0, define subspaces of S'{R) by 

4+ := {/, e^^V e L^}, Li := {/, c^/ G L^}, Li_^+ := {/, a;/ G Lj+j. 

Definition 5.1. Given a constant 77 > and a weight function lo satisfying (1.21), 
define the Banach spaces 

Bi,B2,Xi,X2CV'{R;ClxC,) 

by 

Xi:= Hi X2:=H^^n{d,LlxLl^+). 

with norms 

\\{v,z)\\b, ■■= \\{v,z)\\Hr, 

\\id^v,z)\\B2 ■■= \\id^v,z)\\Hi + \\v\\li + ||e''(')^2:||x,i, 

\\{v,z)\\x, := !|(^,z)||h2, 
II (9x^^,^^)11x2 := II (^a;^^,^) 11^2 + II^^^IIli + ll^e''(')^z||ii. 

In particular, ^2^X1^ Bi, with ||-||bi < |Mlxi < II '11x2 5 ^-iid X2 "-^ B2 "-^ Bi, 
with II • ll^j < II • \\s2 < II • 11x2 5 and the unit ball in Xi is closed in Bi. 

5.1 The perturbation equations 

If solves (1.2) with initial datum U^^^q = U^ + Uq, then the perturbation variable 

U{£,x,t) := f/^(x,t) — U^{x) satisfies 

f dtU - L{e)U = d^Qfie, U,d^U) + Qr{e, U), 
[ U{s,x,0) = U^{x). 

The nonlinear term Qf satisfies (4.2), while Qr satisfies Lemma 4.2. 
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5.2 Coordinatization 



Let <^j_ be the eigenfunctions of L{£) associated with the bifurcation eigenvalues 

'~f{£) lb ir(e), and let (^|. be the corresponding left eigenfunctions. We know from 
Section 3.1.3 that (7 it iT)(e) G C \ C_ U C+, hence decay exponentially at both 
—00 and +00, in particular, if in (1.21) is small enoug h, then € H^. Let U be 
the L^-projection onto span((^j_) parallel to span((pj_)-'-. Decomposing 

U = uu^^% + tfi2^^<y?+ + U2, Uq = oi$R<y9^ + a2Q<p% + b, 

where un^ip^ + Ui2^(p^ and ai^ft^?^ + a2^ip^ belong to span(0^) (so that, in 
particular, uij and aj are real), and coordinatizing as {ui,U2), ui := {uii,ui2) S M^, 
we obtain after a brief calculation that U solves (5.1) if and only if its coordinates 
solve the system 



(5.2) 



where 



dtU2 = (1 - U)L{e)u2 + (1 - n)iV(£, ni,«2), 
%|t=o = a> 
, U2\t=0 = b, 

N{e,ui,U2) := {d^Qf + Qr){e,U' ,U). 



Given Tq > 0, there exist Co > and C > 0, such that, if \a\ + ||6||i^2 < 
Co, the initial value problem (5.2) possesses a unique solution (txi, U2)(a, 6, e) G 
(70([0,ro],M2 xif2) satisfying 

C~Vl - C\\bfH, < \ui{t)\ < C{\a\ + ||6||^2), 
(5-3) ||u2(t)k2 <C(||5|U2 +|a|2), 

\\d(a,b)iui,U2){t)\\c(R2y,H\m) < C. 

(For more details on the initial value problem (5.2) and estimate (5.3), see [TZ2], 
Proposition 4.2.) 

5.3 Poinccire return map 

We express the period map (a, 6, e) b := U2{a,b,e,T) as a discrete dynamical 
system 

(5.4) S = S(£,r)6 + iV(a,6,£,r), 

where 

S(e,T) := e^(i-n)^(^) 
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is the linearized solution operator in v and 

iV(a,6,£,r) := / S{e,T - s){l-Ii)N{e,ui,U2){s)ds 
Jo 

the difference between nonlinear and linear solution operators. 

Evidently, periodic solutions of (5.2) with period T are fixed points of the period 
map (equilibria of (5.4)) or, equivalently, zeros of the displacement map 

A(a, b, s, T) := {S{e, T) - ld)h + N2{a, b, e, T). 

5.4 Lyapunov-Schmidt reduction 

We now carry out a nonstandard Lyapunov-Schmidt reduction following the "inverse 
temporal dynamics" framework of [TZ2], tailored for the situation that Id — S(e, T) 
is not uniformly invertible, or, equivalently, the spectrum of (1 — n)L(e) is not 
bounded away from {jir /T}j^z. In the present situation, (1 — n)L(e) has both a 
1-dimensional kernel (a consequence of (H4), see Section 1.5) and essential spectra 
accumulating at A = 0, and no other purely imaginary spectra, so that Id — S(e, T) 
inherits the same properties; see [TZ2] for further discusssion. 

Our goal, and the central point of the analysis, is to solve A(a, e, T) = 
for & as a function of (a,e,T), eliminating the transverse variable and reducing to 
a standard planar bifurcation problem in the oscillatory variable a. A "forward" 
temporal dynamics technique would be to rewrite A = as a fixed point map 

(5.5) 6 = S(e,r)6 + iV(a,6,e,r), 

then to substitute for T an arbitrarily large integer multiple jT. In the strictly 
stable case ^a{{l — n)L) < < 0, \\^{e,jT)\\iii^Xi) < k for j sufficiently large. 
Noting that N is quadratic in its dependency, we would have therefore contractivity 
of (5.5) with respect to 6, yielding the desired reduction. However, in the absence 
of a spectral gap between a{{l — n)L) and the imaginary axis, ||S(£, jr)||£(j(^^) does 
not decay, and may be always greater than unity; thus, this naive approach does 
not succeed. 

The key idea in [TZ2] is to rewrite A = instead in "backward" form 

(5.6) 6= (Id-S(£,r))-^ iV(a, 6, £,T), 

then show that (Id — S(e, T))"^ is well-defined and bounded on Range TV, thus 

obtaining contractivity by quadratic dependence of N . Since the right inverse 
(Id - S(e, T))"' iV is formally given by ^j^o ^(^li^)-^ ^^^^ amounts to establish- 
ing convergence: a stability/cancellation estimate. Quite similar estimates appear 
in the nonlinear stability theory, where the interaction of linearized evolution S 
and nonlinear source N are likewise crucial for decay. The formulation (5.6) can 
be viewed also as a "by-hand" version of the usual proof of the standard Implicit 
Function Theorem [TZ2]. 
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Lemma 5.2. Under the assumptions of Theorem, 1.14, if the constant rj in Definition 
5.1 satisfies rj < r]o, where rjo was introduced in Corollary 1.6, then 

N : (a, b, e, T) e x Xi x ^ N{a, b, e, T) e X2, 

is quadratic order, and as a map from x ;Bi x to B2 for uniformly 
hounded, with 

\\N\\x,<C{\a\ + \\h\\x,)\ 
(5.7) \\dia,b)N\\c(R^-,B„B,) < C{\a\ + \\b\\H2). 

ll%,r)iV|U(M2,H,) < C{\a\ + IMn^f, 

Proof. We use the variational bounds of [TZ3] (see Propositions 5 and 6, [TZ3]) and 
Lemma 4.2. Note that, in (1.21), only u!~^ G L°° and (1.21) (ii) were used at this 
point. □ 

5.4.1 Pointwise cancellation estimate 

We now develop the key cancellation estimates, adapting the pointwise semigroup 
methods of [ZH, MaZ3, Z2] to the present case. 

Our starting point is the inverse Laplace transform representation (3.55). De- 
forming the contour using analyticity of Q\ across oscillatory eigenvalues A±(e) we 
obtain Q = Q + O, where 

0{x,t-y) := e^+(^)V+(a;)<^+(y)*^ + e^-^^^^ipt{x)<fl{yy'' 

is the sum of the residues of the integrand at X± (the right- and left-eigenfunctions 
<fj_ and (fj^ are defined in Section 5.2). The Green function Q is the kernel of the 
integral operator S{e,t) defined in Section 5.3. Note that, under the assumptions 
of Theorem 1.14, the Evans function associated with (1 — n)L(e) satisfies (1.18), so 
that, by Remark 3.28, Proposition 3.23 applies to Q. 

For I/, vq > 0, let r be the counterclockwise arc of circle dB(0, r) (r as in Propo- 
sition 3.21) connecting —v — ivQ and —ly + ivQ. If ly and j/q > are sufficiently 
small, then F is entirely contained in the resolvent set of (1 — n)L(e), and G can be 
decomposed as -I- Qu, with 



(5.8) 



Qiie, X, t; y):=^.J^ e^'Gxie, x, y) d\ 

Guie, X, t; y) := — P.V. / + / e^'Gxie, x, y) dX 



Let Si and Sn denote the integral operators with respective kernels G\ and Gwi 
so that S = Si -I- Sn, and let := (— £o,£o) x (0) +00); for some £0 > 0. 
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Remark 5.3. The contour T being contained in the resolvent set of L, the elemen- 
tary bound holds: 

\dygx\<Ce-<'^\'-'-y\, Aer, 

for some 9,^ > depending on v. See for instance Proposition 4-4> [MaZ3]. 

Our treatment of the high-frequency term follows [TZ3]: 

Lemma 5.4. Under the assumptions of Theorem 1.14, the sequence of operators 
with kernel '^^=0 ^n(^i '^^) absolutely convergent in C{H^), uniformly in (e, T) G 

n. 

Proof. Starting from the description of the resolvent kernel given in Proposition 
3.21, we find by the same inverse Laplace transform estimates that give terms H 
and TZ in Proposition 3.23, that the high-frequency resolvent kernel Qu, defined in 

(5.8) , may be expressed as 

(5.9) Gil = Ce-^(l^-^l+*) + hT^+stS, 

where C and its space-time derivatives are bounded, 9 > 0, and hTx+stS is a generic 
hyperbohc term; in particular h has the form (3.61) and satisfies (3.64) (i). The 
Lemma follows. □ 

Next we turn to the low-frequency component of Q. Its fluid terms are handled 
as in [TZ3]: 

Lemma 5.5. Under the assumptions of Theorem 1.14, the sequence of operators 
with kernel X^^=o '^-^) converges in C{dxL^ , H^), uniformly with respect to 
(e,T) G n. 

Proof. We argue as in the proof of Proposition 3 of [TZ3]. Let f G L^. By (5.8), 
'l2n=o ^y^^f decomposes into I — Hat, where 

For small u and X E T, (1 - 6^^)-^ = X'^T-^^l + 0(A)). 

The boundary term I is independent of N and is seen to belong to by Remark 

5.3. 

By (3.56)-(3.60), \~^dy£x and X~^dySx have the same form as £\ and Sx. By 
Proposition 3.21, X~^dyTZx behaves like the sum of Tlx ^^nd a pole term of form 
y-i^^e\x-y\^ Hence, by the same Riemann saddle-point estimates used to bound Q 
in Proposition 3.23, we find that 

(5.10) / T^dyGx dX={£ + S + n) (£, ATT), 
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up to a constant (independent of N) term of the form Ce~^l^~^l, where the space- 
time derivatives of C are uniformly bounded. This constant term satisfies the same 
bound as term I. 

In (5.10), £, S, TZ denote generic excited, scattered and residual terms of form 
(3.58), (3.59)-(3.60) and (3.64) (ii). By dominated convergence, 

H^- lim [ £{e,NT)f{y)dy 

exists and is equal to a sum of terms of the form 

(5.11) C{e,T){U'y [ f{y)dy. 

Jr 

Besides, by (3.59)-(3.60) and (3.64), 



(5.12) 



{S + n){e,NT)f{y)dy 



<c(ivr)-/'--i^||/IUi. 



This proves convergence in of the sequence Iljv- □ 

We examine finally the contribution to the series ^„S(e,nT) of the new (not 
present in [TZ3]), reactive terms. 

Lemma 5.6. Under the assumptions of Theorem 1.14, the sequence of operators 
with kernels X^^o absolutely convergent in C{L^+,H^), uniformly 
with respect to {e,T) E il. 

Proof. Let / G L^. By (5.8), Proposition 3.23 and Corollary 3.26, the low-frequency 
kernel Qi satisfies 

(5.13) / e-''y^gi{e,t)e+fiy)dy= f e-^^+ (5 + 7^) (£, t)4/(y) dy, 

Jr jr 

and, by Corollary 3.27, 

/ e-^\s + n){e,t)e+fiy)dy < C(l + ti)e-^i*||/|Ui, 
Jr m 

and the upper bound defines for t = NT an absolutely converging series in i^^. □ 



Prom Lemmas 5.4, 5.5 and 5.6 and the fact that S(e,T) G C{Bi), for ah (e, T) G 
il, we can conclude that, under the assumptions of Theorem 1.14, the operator 
Id — S(£,T) has a right inverse 

(Id-S(£,r))-i :=^S(£,nr) : B2 ^ B^, 

n>0 

that belongs to C{B2,Bi), locally uniformly in (e, T) G 

We will need the following regularity result for the right inverse: 
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Lemma 5.7. Under the assumptions of Theorem 1.14, the operator (Id — S(£:, T)) ^ 
is in {s,T) G fi, with respect to the C{B2,Bi) norm on 82- 

Proof. Note that, by (3.55), dtQ has kernel A^^; in particular, the small A (low- 
frequency) estimates of the proofs of Lemmas 5.5 and 5.6 imply the convergence 
of the sequence X]i^=o ^-^) -^('62, -Bi). The contribution of dTGu{£,nT) is 

handled as in Lemma 5.4, by (5.9) and (3.64) (i) with k = 1. 

Bounds for e-derivatives are handled as in [TZ3], using either the variational 
equation (L — X)dsQ\ = —{dsL)Qx, or the d^Q bounds of Proposition 3.11 from 
[TZ2]. □ 

Note that the e-derivative bound (5.7)(iii) is stated on a proper subspace of Bi, 
namely Xi. In this respect, the following Lemma, asserting boundedness of the right 
inverse on X2 '-^ B2, in C{X2,Xi) norm, is key to the reduction procedure of the 
following Section. (See Remark 5.12.) 

Lemma 5.8. Under the assumptions of Theorem 1.14, (Id — S(e,T))~^ belongs to 
C{X2,Xi), for all {e,T) G n. 

Proof. The convolution bound 



(5.14) 



UJ2 



e-0\-y\f{y)dy 



< Cmin 



L2 



L2 



where C depends on ||a;2e~^l''||£i|-|j;^2, holds by (1.21)(i) and (1.21)(iii). It implies 
that the contributions of Qn, of 1 and of the constant pole term in lljv (see the proofs 
of Lemmas 5.4 and 5.5) are all bounded in C{Xi). 

The scattered and residual terms in II contribute nothing to the limit, by (5.12). 

We use again Corollary 3.27 to handle the contribution of the reactive term. In 
(3.66), there are two terms in the upper bound for {S + 7^)^|. The first term is 
handled by (5.14), and the second by 



(5.15) ( 

noting that 
(5.16) 



-et 



J f e-\^-y\'/^'f{y)dy 
Jr 



< e 



-et 



L2 



^ig-|.-y|2/M* 



L2 



Lh 



uj2e 



-\x~y\yMt, 



L2 < Ce- 



< Ce 



\uj2e 



-Bit 



{tiuj{t)^ 



2 + a;2e' 



-|.|/M| 



L2 



). 



and the upper bound in (5.16) defines for t = NT an absolutely converging series, 
if ^0 in (1-21) is small enough. We used in (5.16) the growth assumption on u. □ 
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Remark 5.9. The above proof shows that, in (1.21), we need in particular 9o < 
where rjQ is the decay rate of the background profile (see Corollary 1.6). In addition, 
we need 60 < mm(0i^,6, -^), where 6 is the rate of decay in the upper bounds of 
Proposition 3.23 and M is the diffusion rate in the upper bound of (3.66) (ii), and 
9q to be small enough so that G H^. 

5.4.2 Reduction 

Corollary 5.10. Under the Assumptions of Theorem 1.14, the equation 
(5.17) A(a,6,£,T)=0 (a,6,e,r) G ^ Xi x 

where A is defined in Section 5.3, is equivalent to 

6 = (Id - S(e, T))-^ N{a, b, e,T)+uj 

for 

w GKer(Id-S(e,r))nXi. 

Proof. A simple consequence of the definition in the above Section of the right 
inverse (Id-S(£,r))-^ For more details, see the proof of Lemma 2.8 in [TZ2]. □ 

Note that, as a consequence of (H4), Section 1.5, the kernel of Id — S(e, T) is of 
dimension one, for all e,T, generated by {U^)'- 

Corollary 5.11. Under the Assumptions of Theorem 1.14, the map 

T{a,b,e,T,a) := {Id - S{e,T)y^ N {a, b,s,T) + a{U^y , 

is bounded from R^ x Xi x R^+i to Xi, and from R^ x Hi x R^+i to Xi, for \a\ 
bounded and \a\ + \\b\\xi + |£| sufficiently small, with 

imU, < C(|a| + IHI^J, 
\\d(a,b)'^\\c{R^xBuBi) < C{\a\ + ||6||xi), 

\\dTT\\B,<C{\a\' + \\bf^J, 
\\da\\B,<C{\a\'' + \\bfx^ + \a\), 
\\daT\\B,<C{\a\' + \\bf^^ + l). 

Proof. Follows from Lemma 5.2, the results of Section 5.4.1, and the above remark 
on the kernel of Id - S (e, T) . □ 

Remark 5.12. Without Lemma 5.8 but with Lemmas 5.4 to 5.7, we could see T 

has a map from Bi to Bi (for fixed a,e,T,a), with quadratic bound \\T\\q^ < 
C||6||g^, and derivative bound WdbTW^f^Q^ < C||6||b^. This would be sufficient to prove 
existence of a fixed point b & Bi as in the following Proposition, but not to prove 
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regularity of the fixed point with respect to e, precisely because the e-derivative bound 
(5.7) (iii) requires more regularity in b, by d^e^^^^^ = e^^^^^d^L (see also Remark 3.12, 
[TZ2]^. Without regularity of the fixed point of T, we could not use the standard 
implicit function theorem in Section 5.4.3. These issues were discussed in detail in 
Section 2.3 o/[TZ2]. 

Proposition 5.13. Under the assumptions of Theorem 1.14, there exists a function 
P{a,e,T,a), bounded from W^'^^ to Xi and from M^+-^ to B\, such that 

A{a,p{a,e,T,a),e,T) = 0, 

\mx, + \\d^e,T)(3\\ciR-^,B,)<C{\a{' + \a\), 

(5.18) \\daP\\ciW^,B,) < C\a\, 

\\daP\\c{R,Bi) - ^' 

for \{a,e,a)\ sufficiently small. Moreover, for |(a,e)|, \\b\\xi sufficiently small, all 
solutions of (5.17) lie on the 1-parameter manifold 

{b = p{a,e,T,a), a G R}. 

Proof. A consequence of the Banach fixed point theorem applied to map T. For 
more details, see the proof of Proposition 2.9 in [TZ2]. □ 

5.4.3 Bifurcation 

The bifurcation analysis is straightforward now that we have reduced to a finite- 
dimensional problem, the only tricky point being to deal with the 1-fold multiplicity 
of solutions (parametrized by a). 
Define to this end 

/9(o, £, T, d) := P{a, e, T, \a\d), 
with a restricted to a ball in M}, noting, by (5.18), that 

ll^llxi, \\d(a,e,T,a)P\\c(R'i,Bi) < C\a\, 

with (3 Lipshitz in (a, e, T, d) and away from o = 0. 
Solutions {ui,U2) of (5.2) originating at 

(5.19) {a,b) = {aJ{a,e,T,a)), 
by (5.3), remain for < t < T in a cone 

C := {{ui,U2) : |tt2| < C'll'Uil}, 
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Ci > 0. Indeed, (5.3) implies the bound 

\\u2mx, < C\ui{t)\, for ||6|Ui < Ci|a|, 

for Ci > small enough, for all < t < T. 

Likewise, any periodic solution of (5.2) originating in C, since it necessarily 
satisfyies A = 0, must originate from data (a, 6) of the form (5.19). 

Defining b = (3{a, e, T, a), and recalling invariance of C under flow (5.2), we may 
view v{t) as a multiple 

(5.20) U2{x,t) = c{a,£,T,a,x,t)ui{t) 

of ui(t), where c is bounded, Lipschitz in all arguments, and away from a = 0. 
Substituting into (5.2) (i), we obtain a planar ODE 

= ( ^(J)) "1 + M{ui,e, T, t, a, a) 

in approximate Hopf normal form, with nonlinearity M := UN now nonautonomous 
and depending on the additional parameters {T,a,a), but, by (4.2) and (4.4), still 
satisfying the key bounds 

(5.21) |M|, |5,,T,dM| < C|ni|2; |5„,^M| < C|ui| 

along with planar bifurcation criterion (1.22). Prom (5.21), we find that M is in 

all arguments, also at a = 0. By standard arguments (see, e.g., [HK, TZl]), we thus 
obtain a classical Hopf bifurcation in the variable ui with regularity C^, yielding 
existence and uniqueness up to time-translates of a 1-parameter family of solutions 
originating in C, indexed by r and S with r := ai and (without loss of generality) 
a2 = 0. Bound (1.23) is a consequence of (5.3) (i) and (5.20). 

Finally, in order to establish uniqueness up to spatial translates, we observe, 
first, that, by dimensional considerations, the one-parameter family constructed 
must agree with the one-parameter family of spatial translates, and second, we 
argue as in [TZ2] that any periodic solution has a spatial translate originating in 
C, yielding uniqueness up to translation among all solutions and not only those 
originating in C; see Proposition 2.20 and Corollary 2.21 of [TZ2] for further details. 

6 Nonlinear instability: Proof of Theorem 1.15 

We describe a nonlinear instability result in general setting. Consider 
(6.1) dtU = LU + d^N{U) + R{U), 

well-posed in iJ*, where 

L = d^{Bd^U) + d^U) + GU, 
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and \N{U)\, \RiU)\ < (7|?7p for \U\ < C. Suppose that L has a conjugate pair of 
simple unstable eigenvalues A± = 7 it fr, 7 > 0, and the rest of the spectrum is 
neutrally stable, without loss of generality e^^~^^^^ < Ct, where 11 is the projection 
onto the eigenspace associated with X±. 

Coordinatizing similarly as in Section 5 by 

U{x,t) = un(fi{x) + ^12^^2(2;) + U2ix,t), 

where (fj = 0(e~^l^l) are eigenfunctions of L, denote r{t) := \ui\{t). Then, so long 
as |J7|_ffs < Ce, we have existence (by variation of constants, standard continuation) 
of solutions of (6.1) in H^, with estimates 

.g2^ r' = jr + 0{e)\U\, 

^ ' ^ «2 = ii-n)Lu2 + o{€)\u\ 

in L\ 

We shall argue by contradiction. That is, using (6.2), we shall show, for C > 
fixed, e > sufficiently small, and \u2{0)\hi < Cr(0), that eventually r{t) > e, no 
matter how small is r(0), or equivalently \U\h^{0)- This, of course entails nonlinear 
instability. 

Define a{x,t) := U2{x,t)/r{t). Then, 

f ru'2 — U2r' _ u'2 U2 r' 

yielding after some rearrangement the equation 

(6.3) a' = ((1 - U)L - 7)a + 0(e(e-^l^l + \a\ + |ap)) . 

From (6.3) and standard variation of constants/contraction mapping argument, we 
find that |a(t)|/^i remains less than or equal to C|a(to)|iji for t — to small. 

By variation of constants and the semigroup bound |e'''^^~^^^""'''^*|^i_,.^i < Ce~'^ 
(note: 7 is scalar so commutes with (1 — n)L), we obtain 

<5(t)<C(|a(0)|^i+e(l + 5(i)2)), 

for 5{t) := supo<,-<f |a(r)|j:/i. So long as S remains less than or equal to unity and 
Ce < ^, this yields 

S{t)<2C{\a{0)\Hi+e), 

and thus 

,5(t)<2C|a(0)|^i+^. 
Substituting into the radial equation, we obtain 

r' > (7 - (1 + 6)e)r, 
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yielding exponential growth for e sufficiently small. In particular, r > Ce for some 
time, and thus |?7|j/i > e, a contradiction. We may conclude, therefore, that \U\x,2 
eventually grows larger than any e, no matter how small the initial size r(0), and 
thus we may conclude instability of the trivial solution U = 0. Taking now (6.1) 
to be the perturbation equations about a strong detonation profile, we obtain the 
result of nonlinear instability of the background profile U. 

Remark 6.1. In the easier case of a single, real eigenvalue, the scalar, w equation, 
would play the role of the radial equation here. This case is subsumed in our analysis 
as well. 
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